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We present the Usadel theory describing the superconducting proximity effect in heterostructures 
with a half-metallic layer. It is shown that the full spin polarization inside the half-metals gives 
rise to the giant triplet spin-valve effect in superconductor (S) - ferromagnet (F) - half-metal (HM) 
trilayers as well as to the <po-junction formation in the S/F/HM/F/S systems. In addition, we 
consider the exactly solvable model of the S/F/HM trilayers of atomic thickness and demonstrate 
that it reproduces the main features of the spin-valve effect found within the Usadel approach. 

Our results are shown to be in a qualitative agreement with the recent experimental data on the 
spin-valve effect in MoGe/Cu/Ni/Cr02 hybrids [A. Singh et al, Phys. Rev. X 5, 021019 (2015)]. 

PACS numbers: 74.62.-c, 74.78.Fk, 74.45.+C, 72.25.-b 


I. INTRODUCTION 

Spin-polarized superconducting states attracts grow¬ 
ing interest since they are expected to provide power¬ 
ful mechanisms for controlling the current and magne¬ 
tization in the devices of superconducting spintronicsii 
Although polarized states are not supported in conven¬ 
tional s-wave superconductors, they can emerge in ar¬ 
tificial heterostructures consisting of a superconductor 
(S) and several ferromagnetic (F) layers with different 
orientations of magnetic moments.- 2-5 The non-collinear 
exchange field in the F-layers destroys the spin-singlet 
structure of Cooper pairs penetrating from the super¬ 
conductor. This results in the appearance of spin-triplet 
superconducting correlations with all possible spin pro¬ 
jections S , both polarized states with S = ±1 and non¬ 
polarized one with S = 0. 

The correlations with Ail have two distinctive fea¬ 
tures which give an insight into their experimental ob¬ 
servation and practical utilization^^ First, such corre¬ 
lations are unsensitive to the exchange field parallel to 
the spin quantization axis. As a result, they become 
long-range: in diffusive systems their decay length inside 
the ferromagnet is comparable with the one in normal 
metal while the correlations with S = 0 decay at much 
shorter distances from the superconductor. Second, these 
long-range triplet correlations (LRTC) appear only if the 
ferromagnet has a non-collinear distribution of magneti¬ 
zation. Thus, to control the amplitude of the LRTC in 
the system one can use ferromagnetic bilayer with tun¬ 
able mutual orientation of magnetic moments while the 
effects coming from the correlations with S = 0 can be 
damped by increasing the F-layer thickness. 

Experimental observation of the LRTC is mostly based 
on probing the long-range Josephson current in the 
S/F'/F/F'/S junctions^ - — If the thickness of the cen¬ 
tral F layer is believed to exceed the decay length of the 
short-range non-polarized correlations the observation of 
non-zero critical current can be attributed to the pres¬ 
ence of the LRTCJ^ Other observations of the LRTC are 


based on the so-called triplet spin-valve effect in S/Fi/F 2 
and Fi/S/F 2 systems revealing in non-monotonic depen¬ 
dence of the S-layer critical temperature T c on the an¬ 
gle 0 between the magnetic moments in the Fi and F 2 
layers^ - — 7 The LRTC open an additional channel for the 
“leakage” of the Cooper pairs from the superconductor. 
As a result, T c can have the minimum at 9 ^ 0, 7 T i 18 i 19 
However this effect in T c (0) is typically washed out by 
the monotonically increasing contribution from the cor¬ 
relations with S = 0, which shifts the minimum of T c 
from 9 = 7t/2. 

During the past few years the focus in the studies of 
LRTC is moving towards the heterostructures contain¬ 
ing half-metallic (HM) layers (e.g. CrO 2 ) ] 20 i 21 A recent 
progress in fabrication of such structures has resulted in 
several breakthrough experiments manifesting the long- 
range Josephson current through the layer of Cr0 2 22 ’ 23 
and triplet spin-valve effect in MoGe/Cu/Ni/Cr02 
structures. 24 - The importance of these experiments is con¬ 
nected with the fact that in half-metals the energy bands 
for electrons with spin up and down are separated at a 
distances comparable with the Fermi energy. As a result, 
only spin-polarized correlations with S = +1 can pene¬ 
trate into the HM-layer while all other correlations should 
vanish at its boundary and cannot influence the Joseph¬ 
son current in S/F/HM/F/S junctions or the triplet spin- 
valve effect in S/F/HM systems. Thus, half-metals pro¬ 
vide a unique possibility to probe the phenomena caused 
by LRTC independently from other effects. 

However up to now there is no convenient and com¬ 
monly accepted theoretical model describing the su¬ 
perconducting proximity effect with half-metals. Most 
common approaches for the treatment of the proxim¬ 
ity effect in multilayered structures are based on the 
quasiclassical approximation^^ This approximation be¬ 
comes broken near the interfaces which implies using 
some sort of boundary conditions matching the quasi¬ 
classical Green functions in different layers. In con¬ 
trast with the S/F interfaces where near the critical 
temperature one may use linear Kupriyanov-Lukichev 
boundary conditions for the anomalous Green function^ 
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the interfaces with half-metals require more sophisti¬ 
cated boundary conditions since the number of the Green 
function components at the opposite sides of the inter¬ 
face is different due to the large energy separation of 
the spin-up and spin-down bands. There are a lot of 
papers where the authors made attempts to overcome 
this problem using different versions of the scattering 
matrix approach. The resulting boundary conditions 
were extensively used for the description of half-metals 
within the Blonder-Tinkham-Klapwijk)2£r— Bogoliubov- 
de-GenneS)22r— Eilenberger^r— and Usadel^r— for¬ 
malisms. Despite the fact that these models provide a 
number of generic qualitative predictions, most of the re¬ 
sults strongly depend on the microscopical mechanisms 
of the singlet-triplet conversion and on the concrete form 
of the scattering matrices ) 35 ' 42 ’ 43 which are not available 
form the experimental data. An alternative phenomeno¬ 
logical approach based on the circuit theory^ does not 
contain any information about the particular geometry 
of the system and, thus, can hardly be applied for the 
quantitative description of real heterostructures. 

In the present paper we propose the phenomenologi¬ 
cal model of the superconducting proximity effect with 
half-metals based on the Usadel equation in the diffu¬ 
sive limit. Our model is based on the following three 
key assumptions, (i) The impurities are non-magnetic 
and do not cause spin-flip processes, which allows to in¬ 
troduce the anomalous Green function inside the half¬ 
metal. The Green function component with S = +1 sat¬ 
isfies the Usadel equation equivalent to the one in normal 
metal while all other components are zero, (ii) The are 
no barriers at the boundaries of the HM layers which re¬ 
sults in the continuity of the component with S = +1. 
(iii) The components with S = 0 and S = — 1 cannot 
penetrate the HM layer and vanish at its outer bound¬ 
aries. Our simple model is shown to explain all main 
features of the recently observed triplet spin-valve effect 
in MoGe/Cu/Ni/Cr0 2 structures^ and predict several 
unusual phenomena manifesting the differences between 
the influence of weak and strong ferromagnets on the 
proximity effect. 

To verify our main conclusions we also considered the 
situation when the layers of the S/F/HM spin-valve have 
atomic thickness and are separated by the tunnel bar¬ 
riers. The advantage of such microscopical model is the 
possibility to find the exact solution of the Gor’kov equa¬ 
tions without any prior assumptions about the profiles of 
the Green functions. Using this model we obtained the 
analytical dependencies T c (6) which appear to reproduce 
all main features found within the Usadel approach. 

The paper is organized as follows. In Sec. [IT] we intro¬ 
duce our model. In Sec. m we apply it for the descrip¬ 
tion of the triplet spin-valve effect in S/F/HM systems 
and compare our results with the experimental observa¬ 
tions of Ref. HH In Sec. [TV] we study the anomalous 
Josephson effect in dirty S/F/HM/F/S structures with 
non-coplanar magnetic moments and demonstrate that 
such systems support the states with the spontaneous 



FIG. 1: (Color online) The sketch of the S/F/HM spin valve. 
The exchange field in the F-layer makes the angle 9 with the 
spin quantization axis in the half-metal. 

phase difference, previously predicted for the ballistic 
limit.— In Sec. |V] we consider the spin-valve effect for 
the S/F/HM systems of atomic thickness and compare 
the results with the conclusions of Sec. urn In Sec. IVIl we 
summarize our results and discuss their possible applica¬ 
tions. 


II. MODEL 

Let us consider a multilayered structure consisting of 
superconductors, ferromagnets and half-metals with the 
interfaces perpendicular to the axaxis. Two examples 
of such structures are shown in Fig. |T] and Fig. [4j We 
assume that the systems is in the diffusive limit and the 
temperature T is close to the critical temperature T c of 
the superconducting transition. In this case outside the 
half-metallic layers the superconducting properties of the 
system can be described in terms of the linearized Usadel 
equation^ 

+ + A = 0, (1) 

where the quasiclassical Green function 

/= (/U /It) =V' + tt& ') i& v (2) 

is the 2x2 matrix in the spin space, A = A ia y is 
the superconducting pairing potential inside the S lay¬ 
ers, = 7rT(2n +1) > 0 are the Matsubara frequencies, 
h is the exchange field in the ferromagnets, and D is the 
diffusion constant. 

To describe the superconducting correlations inside the 
half-metallic layers we assume that (i) electron scatter¬ 
ing on the impurities does not cause spin-flips and (ii) 
the barriers are not spin-active. In this case we may in¬ 
troduce the Usadel Green function with only one nonzero 
component meaning that the spin polarization in the 
half-metal is directed along the 2 -axis. 

To match the solutions of the Eq. © in different layers 
one should put the boundary conditions at each interface. 
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We assume all interfaces between the layers to be trans¬ 
parent for electrons. In this case at the interfaces which 
separate two non-halfmetallic material the boundary con¬ 
ditions come to the continuity of / and the combination 
ad x f (cr is the Drude conductivity of the corresponding 
layer) Similarly, for the interfaces with the half-metal 
we demand only the continuity of fff and crd x f- 1 -|- while 
all other components fn, and /m|- should vanish at 
the HM-layer boundaries. At the outer boundaries of the 
heterostructure we demand d x f = 0. 

The physical meaning of the boundary conditions at 
the interface with half-metal is very clear: the interface 
with the HM layer plays the role of spin filter, which is 
absolutely transparent for the correlations with S = +1 
and opaque for all other correlations. 


III. TRIPLET SPIN-VALVE EFFECT IN 
S/F/HM SYSTEMS 

The geometry of the S/F/HM system under consider¬ 
ation is shown schematically in Fig. [1] We choose the 
origin of the rr-axis in a way that the superconductor is 
at — d s < x < 0, the ferromagnet is at 0 < x < df , and 
the HM layer occupies the region df < x < df + dh- The 
exchange field h in the F-layer is assumed to be rotated 
on the angle 6 in the xz-plane and have two components: 
h z = h cos 6 and h x = hsmd. The spin quantization axis 
in the HM-layer coincides with the z one. 

To demonstrate the key difference between the 
S/F/HM system and classical S/F 2 /F 2 spin valves it is 
instructive to write the explicit form of Eq. CD inside the 
F-layer for all components f s and f f : 

(Df/2)dlf s = uj n f s + ih cos 9 f tz + ihsm0f tx , 

(Df/2)dlf tz = ix n ftz +ihcos0f s , 

( 3 ) 

(Df/2)d x f tx = wjd + ihsinOfs, 

{Df/2)d x ft v = IX n fty , 

where Df is the diffusion constants in the ferromagnet. 
Clearly the spin-singlet component f s of the anomalous 
function coming from the superconductor induces the 
triplet components ft z and ft x while the equation for 
ft y remains independent. In the S/F 2 /F 2 structures the 
boundary conditions do not mix different components of 
the function / and, thus, the absence of the source in 
equation for f ty immediately leads to ft y {x) = 0 in the 
whole heterostructure. 

However, in the S/F/HM systems the situation is com¬ 
pletely different. The penetration of the non-zero com¬ 
ponent fff = —ftx + 7 fty into the HM-layer together with 
the vanishing of the component /jq = ft x +ift y is possible 
only if inside the half-metal 


Consequently, this results in the appearance of f ty also 
in the S and F layers. 

For the further analysis it is convenient to exclude f ty 
from CD by substituting the solution of the equation for 
fty into the boundary conditions. As a result, all infor¬ 
mation about the component f ty becomes included into 
the effective boundary condition for f tx ■ which reads (see 
Appendix lAl) 


d x fu 


= -9/ r 


x=df 


with 


r — 2 fih + 


A+ +A if 
1 + fi 8 nf' 


In this expression 


H = 



tanh^dj), 


( 5 ) 

( 6 ) 


( 7 ) 


where the index j £ {s, /, h} corresponds to the S, F and 
HM layers respectively, qj = -^2ui n /Dj, Dj and cTj are 
the diffusion constant and the normal conductivity in the 
j-th layer. 

Note that in the case when the diffusion constants and 
conductivities of all layers are equal to each other the 
expression for T takes the form 


T = 2tanh (q h d h ) + tanh [qf(d s + df)] . (8) 

Clearly, this expression reflects the fact that at the F/HM 
interface the component f tx of the Green function induces 
two components ( f tx and f ty ) in the half-metal and also 
the component f ty in the S/F bilayer. 

Now let us analyze the dependence of the S-layer crit¬ 
ical temperature on the angle 9. To simplify the cal¬ 
culations we assume that the thickness of the S-layer is 
much smaller than the superconducting coherence length 
£ s o = \]D s /A-kT c o, which enables neglecting the spatial 
variation of the pairing potential A across the supercon¬ 
ducting film (here T c $ is the critical temperature of the 
isolated superconductor). For convenience we choose A 
to be real. Then the dependence T c (0) is defined by the 
self-consistency equation 


A In 


To (9) 

T c o 


+ E 


n —0 


A 

n+ 1/2 


27 tT c ( 9 ) fs ( 9) 


0. (9) 


The equation (ID with the boundary conditions de¬ 
scribed above allows us to obtain the analytical expres¬ 
sions for the anomalous Green function in each layer (the 
details of the calculations are presented in Appendix [Bj. 
Inside the superconductor the component f s has the form 

A \ A HsW 2 sin 2 9 

Era + "TV 1 / Uj n Q sin 2 9 + T + (T -2fi h ) cos 2 9 ’ 

( 10 ) 


fty = iftx ^ 0. 


( 4 ) 
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FIG. 2: (Color online) The dependencies of the critical tem¬ 
perature on the angle 9 for different thicknesses df of the 
ferromagnetic layer. The system parameters are: d a = 0.5£ a o, 
h/AnT c o = 2, dh = \ZD h /4nT c0 , D f /D s = af/a 3 = 1CT 3 , 
and Dh/Df = oh/af = 1CT 3 . 


FIG. 3: (Color online) The dependencies of the critical tem¬ 
perature on the thickness of the ferromagnetic layer for 9 = 0 
(blue curve) and 9 = 7r/2 (red curve). The system parameters 
are the same as in Fig. [2] 


where q = (uj n + ih) /Df , p = p s qf, 


IV. ANOMALOUS JOSEPHSON EFFECT IN 
S/F/HM/F/S SYSTEMS 


W = 3 


_ 1 _\ 

qcosh(qdf) + psmh(qdf) J 


( 11 ) 


J q p + qta,nh(qdf) \ 
W q + ptanh(qd f ) j ’ 


( 12 ) 


and the pair-breaking parameter 


T n 1 = ^ 7 TT<lCOth(qdf) (13) 

(j s Zd s 


is the same as for the S/F/S junction with the F-layer 
thickness 2 df in the 7r-state (see, e.g., Ref. 0 ). 

The expression m allows to analyze the main features 
of the critical temperature behavior. First, in contrast 
with the case of S/F 1 /F 2 spin valve^S the dependence 
T c {9) is always symmetric, i.e. T c {n — 9) = T C (Q). Sec¬ 
ond, f s (sin 2 9 ) is a monotonically decreasing function 
and, thus, the minimum of the critical temperature cor¬ 
responds to 9 = 7t/2. Note that these two features were 
clearly observed in recent experiments with spin valves 
containing the half-metallic CrC >2 layers 24 

The typical dependencies T C (Q) are shown in Fig. [2] 
One sees that S/F/HM system with df ~ £/ reveal gi¬ 
ant triplet spin-valve effect originating due to the LRTC: 
for the chosen parameters the increase of 9 results in the 
damping of T c from 0.4T c o to zero. In Fig. [3] we also plot 
the dependencies of T c on the F-layer thickness df for 
9 = 0 and 9 = 7t/2. For df -C £/ the superconductiv¬ 
ity is fully suppressed by the proximity with half-metal. 
When increasing df above a certain threshold, which de¬ 
pends on 9 , the critical temperature becomes non-zero 
and grows rapidly. Interestingly, the difference in the 
thresholds for 9 = 0 and 9 = n/2 makes it possible to si¬ 
multaneously reach the absolute maximum of T c for 9 = 0 
and have T c (n/2) = 0. 


In this section we analyze the behavior of the Joseph- 
son current through the Fi/HM/F 2 -trilayer (see Fig. [4]). 
In particular, we show that the peculiar mixing of differ¬ 
ent spin-triplet / components at F/HM interfaces results 
in the tpo-junction formation provided the magnetic mo¬ 
ments in the F-layers and the spin quantization axis in 
the HM-layer are the non-coplanar vectors. Previously 
this effect was noted within the circuit-theory for the 
junctions consisting of two ferromagnetic superconduc¬ 
tors separated by the half-metal.— The model of the fer¬ 
romagnetic superconductors^ allows to simplify the cal¬ 
culations but can hardly be applied for the real systems in 
which the ferromagnetic order strongly suppresses super¬ 
conductivity. Here we consider a more realistic situation 
when the regions with superconducting and ferromag¬ 
netic orders are separated in space. 

To make the physical origin of the effects under con¬ 
sideration more transparent we restrict ourselves to the 
simplest case when the exchange field in the Fj-layer is 
directed along the x-axis while the magnetic moment of 
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FIG. 4: (Color online) The sketch of the S/F 1 /HM/F 2 /S 
Josephson junction. The exchange field vectors in the fer- 
romagnets and half-metal are perpendicular to each other. 
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the F 2 -layer has only y-component. For simplicity we 
consider equal magnitude h of the exchange field in the 
Fi and Fi-layers. Similar to the previous sections we as¬ 
sume that in the half-metal the spins are quantized in the 
^-direction. We choose the £-axis perpendicular to the 
interfaces so that the Fi, HM and F 2 -layers occupies the 
regions — {df + dh/2) < x < —dh/2 , —dh/2 < x < dh/2 
and dh/2 < x < df + dh/2 respectively. 

Before proceeding with the calculation of the Joseph- 
son current let us briefly point out the main difference of 
the described system from the Si/F 1 /F'/F 2 /S 2 junctions 
with the same magnetic configuration extensively studied 
before. In the latter system if the magnetic moments in 
the three ferromagnets are perpendicular to each other 
and the thickness of th e centr al layer strongly exceeds 
the coherence length yjD-i/hi the Josephson current is 
negligibly small. Indeed, the Fi layer produces the com¬ 
ponent ftx of the anomalous Green function which be¬ 
comes long-range in the F / layer and reaches the ferro- 
magnet F 2 . However, this component do not contribute 
to the Josephson current since the magnetic moment in 
the F 2 directed along the y-axis cannot convert this ftx 
component into the singlet f s one. 

The situation becomes completely different if instead 
of the ferromagnet F' one has half-metallic layer. In this 
case at the Fi/HM interface both f tx and ft v compo¬ 
nents are produced. As a result, in the F 2 -layer the ft y 
component can be effectively converted into the f s one 
and, thus, produce the non-vanishing Josephson current. 

The current-phase relation j((p) of the Josephson junc¬ 
tion is defined by the sum 

j (<p) = SneDNT £ 3 {f* s d x f s - t/d x f t ) , (14) 

OJn> 0 

where e is the charge of electron, N is the electronic 
density of states, D is the diffusion coefficient. Since 
the current does not depend on the position across the 
junction the anomalous Green function can be taken in 
arbitrary point. Practically it is convenient to choose this 
point, e.g., at the S/Fi interface (at x = —df — dh/2) 
where only the singlet component f s is non-zero. 

To simplify the further calculations we assume that the 
normal conductivity of the superconducting electrodes 
strongly exceed the ones in the ferromagnets, so that at 
both S/F interfaces the rigid boundary conditions are 
fulfilled: f s = (A /ui n )e ± ' lv ^ 2 (the signs + and — in the 
phase factor correspond to the right and left S-layers re¬ 
spectively) while f t = 0. Also we assume that the thick¬ 
ness of the HM-layer is much less than the coherence 
length yjDh /47tT c o which allows to neglect the spatial 
variations of the function /-pj- across this layer. Solving 
the Usadel equation with the boundary conditions dis¬ 
cussed above we obtain the analytical expression for the 
derivative d x f s at x = — df — dh/2 (see Appendix [Cl) and 
the resulting current-phase relation: 

j (tp) = AneDfNTsin (ip + ^ (15) 

V 2 / w „>o w - 



FIG. 5: (Color online) S/F/HM spin valve of the atomic thick¬ 
ness. The angle between the exchange field in the ferromagnet 
and the spin quantization axis in the half-metal is denoted as 
9. The coupling between the layers is described by the trans¬ 
fer integral t. 

where 

Jn = _1_32 [_ q _' 

n qfCoih(qfdf) + $t[qcoth(qdf)] [sinh(gdj) 

(16) 

and the definition of q and < 7 / are the same as in Sec. IIIII 
From Eq. (fl6l) one can see that the non-coplanarity of 
the magnetic moments in the magnetic layers results in 
the appearance of the spontaneous Josephson phase dif¬ 
ference (so-called </?crjunction). In contrast with the ordi¬ 
nary S/F/S systems where the formation of </?o-junction 
requires strong spin-orbit, coupling^ - — here the sponta¬ 
neous phase arises due to spin polarization in the HM- 
layer. Note that in case of arbitrary mutual orientation 
of the magnetic moments in the F-layers one can expect 
the spontaneous Josephson phase <p 0 to be equal to the 
angle between the projections of these magnetic moments 
to the xy -plane (see Ref. [3) . 


V. S/F/HM SPIN VALVE OF THE ATOMIC 
THICKNESS 

In this section we study the spin-valve effect in the 
S/F/HM trilayers of atomic thickness. Experimentally 
such kind of systems can be realized, e.g., on the basis 
of the RuSi' 2 GdCu 2 C >8 and Lao. 7 Cao.sMn 03 compounds. 
The first one is the ferromagnetic superconductor with 
the alternating S-layers of CuC >2 and the magnetically 
ordered Ru 02 layers. The second compound is a half- 
metal^ which is recently shown to have a strong influence 
on the properties of the superconducting systems.— 
Here we use the microscopical Gor’kov formalism to 
calculate how the critical temperature of the S/F/HM 
trilayer depends on the angle 9 between the exchange 
field h = h cos 9z + h sin 0 x in the ferromagnet and the 
spin quantization axis z in the HM-layer (see Fig. GO). 
We assume that the neighboring layers are coupled by 
the electron tunneling describing by the transfer integral 
t. Previously this approach was successfully applied for 
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the description of the spin-valve effect in F/S/F struc¬ 
tures (see [50M53l| and references therein). The main ad¬ 
vantage of the atomic layers model is the possibility to 
obtain the exact solutions for the Green functions and 
analyze their properties. Below we obtain such solutions 
for the case of the S/F/HM structures and demonstrate 
that they reproduce all main features of the spin-valve ef¬ 
fect described within the phenomenological Usadel model 
in Sec. IIIII 

Let us denote the two-component electronic operators 
in the S, F and HM layers as ip, ip and 17 respectively. 
For simplicity we assume that the quasiparticle motion in 
the plane of the S and F-layers is described by the same 
energy spectrum £(p). At the same time, in the half¬ 
metal the energy is strongly spin dependant: we assume 
that for the spin-up quasiparticles £ 4 - = £(p) while for 
the spin-down ones £ 4 . = + 00 , which implies 774 . = 0 . 

The Hamiltonian of the system under consideration 
has the form^S 


H = H 0 + Hs + H t 


where 


(17) 


Ho = + ip+C(e)i/j a +r)+Pr] a , (18) 

p;a=l,2 


Hs = ( A Vp,2V?-p,i + A +p,i+± P , 2 ) > (19) 

P 


Ht= ^2 t (vi^a + i>£<Pa + 'tp£'n a + r]+'llj a ) ■ ( 2 °) 

p;a=l,2 

In Eq. (1181) we introduces two matrices describing the 
effect of Zeeman coupling in the F-layer and the spin 
polarization in the half-metal^ 


r(n\ - ( £ ~ hcose ~h sin 0 \ p_(€ 0 \ 

' y —h sin# £ + hcos9 J ’ y 0 +00 J ' 

( 2 !) 

Note that in contrast with the model of Ref. [441 in our 
system the superconducting and ferromagnetic regions 
are separated in space which makes it possible to consider 
the exchange field h of arbitrary magnitude. 

The critical temperature 7/ of the S-layer is defined 
by the linear expansion of the anomalous Green func¬ 
tion F+p = (r T over the gap potential |A|, 

where T t denotes the time-ordered product for the imag¬ 
inary time r. Writing and solving the systems of Gor’kov 
equations for the S/F/HM system we find (the details of 
calculations are presented in Appendix [Dj: 


F+ 

A* 


j(*w + £)i 


t 2 (iu + C) 


t 2 (ico + P) 1 


-1 


To simplify the further calculations we assume the tun¬ 
neling constant t to be small and perform the power ex¬ 
pansion of Eq. (E^l) over t 2 . To obtain the non-trivial de¬ 
pendence 7/(0) we should keep the terms up to t 6 . Also it 
is convenient to represent the self-consistency equation^ 
in the form 


+ OO 


TJff\ = TJO 1 


p+ (n\ — f+ <TP 


co n >0 


(23) 

where 1/(0) = T c0 [l — 0(t 2 )] is the critical temperature 
at 0 = 0 [for t <C 7 /0 one has |7/(0) — T c0 | -C T c0 ] and the 
sum is taken over the discrete set of positive Matsubara 
frequencies u n = 7rT/o(2n +1). 

Substituting the expression for Fp 2 into Eq. (l23l) we 
found: 


+ OO 


7/(0) = 7/(0)+ 


uj n >Q 


(io 


47' 2 0 t 6 /i 2 sin 2 0d£ 

0(iw + 0 3 [(iw + 0 2 


Taking the integral over £ we get: 

7/(0) = T c ( 0) - 


7 T T 2 0 t 6 h 2 sin 2 0 


2 -—' u; 3 (4u ; 2 
<j„>o v 


h 2 ) 3 ’ 


-h 2 ] 3 ’ 

(24) 

(25) 


First, one can clearly see that the deviation of the crit¬ 
ical temperature from T/(0) is proportional to sin 2 9 and, 
thus, T c (n — 9) = T c {9) in the full accordance with the 
conclusion of Sec. mu 

Second, Eq. (l25l) shows that the magnitude of the 
spin-valve effect which can be characterized by the value 
ST C = 7/(0) — 7/(0) has non-monotonic dependence on h 
with the maximum §T™ ax oc t 6 /T 3 0 at h ~ T c q. Indeed, 
for h < 7 /0 the exchange field weakly affect the system 
properties, and 


_ 511C(9 )t 6 h 2 t 6 f h \ 2 

2 1 WJ 0 ^TfVTW 


(26) 


In the opposite limit when h 7 /0 the strong Zeeman 
splitting of the energy bands inside the F-layer effectively 
damps the tunneling constant between the layers, and as 
a result the spin-valve effect is also weak: 


<57/ = 


7C(3 )t 6 ^ t 6 

8tt 2 7/o/i 4 K T [ 0 



(27) 


In terms of the previously discussed Usadel theory the 
latter relation simply reflects the fact that for h 7 /0 
the coherence length £/ = yjDf/h in the F-layer be¬ 
comes much smaller than its thickness, and supercon¬ 
ducting correlations do not reach half-metal. 


VI. CONCLUSION 


xJ 


(iui 


Oi ~t 2 



C) — t 2 (ico 



-1 

To sum up, we proposed the phenomenological Usadel 
theory of the superconducting proximity effect in multi- 
(22) layered systems with a half-metallic layer. It is shown 
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that the boundary between ferromagnet and half-metal 
serves as a source of additional triplet component of the 
anomalous Green function which is perpendicular to both 
exchange field h in the F-layer and the spin quantization 
axis in half-metal z. For the S/F/HM trilayes we ana¬ 
lyzed the dependence of the critical temperature T c on 
the angle 9 between h and z and found that the discov¬ 
ered triplet component strongly enhances the spin valve 
effect compared to the traditional S/F/F structures: in¬ 
creasing the angle 9 one can damp T c from the value com¬ 
parable to the critical temperature of the isolated super¬ 
conductor down to zero. Note that the described giant 
damping of T c appears only due to the long-range triplet 
correlations (LRTC) since the short range ones do not 
penetrate the HM-layer and, thus, are not sensitive to 9. 
In addition, we showed that the full spin polarization in 
the HM-layer requires the symmetry T c (tt — 9) = T c {9), 
which was clearly observed in recent experiments with 
the MoGe/Cu/Ni/CrC >2 spin valves.— 

To verify our main conclusions about the peculiarities 
of the spin-valve effect in the S/F/HM structures we con¬ 
sidered the case when the layers of such system have the 
atomic thickness. For this case we obtained the exact an¬ 
alytical solution of the Gor’kov equations and calculated 
the dependencies T c {9). We found that if the tunneling 
rate between the layers is small the deviation of T c {9) 
from the critical temperature at 9 = 0 is proportional 
to sin 2 9 which reproduces the symmetry relation found 
within the Usadel formalism. 

Also we demonstrated that the new “perpendicu¬ 
lar” triplet component of the anomalous Green func¬ 
tion dramatically modifies the current-phase relation of 
the S/F/HM/F/S Josephson junctions provided the ex¬ 
change field vectors in the F-layers and the spin quanti¬ 
zation axis z in half-metal are non-coplanar. First, we 
found that such systems support the ^o-junction for¬ 
mation. This result is non-trivial since in the usual 
S/F/S and S/N/S structures the appearance of the spon¬ 
taneous Josephson phase difference requires strong spin- 
orbit coupling.— In contrast, here such phase emerges 
only due to the spin selectivity of the half-metal. Second, 
the critical current of the S/F/HM/F/S structures does 
not vanish when the exchange field vectors and z are 
perpendicular to each other. This result directly origi¬ 
nates from the presence of the additional Green function 
component. It is exactly this component, which makes it 
possible not only to generate the long-range triplet corre¬ 
lations near the S-lead but also to convert them back into 
the singlet ones near the opposite lead. Note that previ¬ 
ously these effects were discussed for the system of two 
ferromagnetic superconductors and the half-metal sepa¬ 
rated with the tunnel barriers— However the model of 
ferromagnetic superconductor is valid only for extremely 
weak exchange field values which make it inapplicable to 
the real heterostructures. 

Finally, we would like to mention that the appearance 
of the additional triplet component of the anomalous 
Green function should strongly influence the local density 


of states (LDOS) and the electromagnetic response of all 
considered heterostructures. Indeed, spin-triplet correla¬ 
tions make positive contribution into the LDOS 44 ' 56- — 
and, as a result, LDOS in the superconducting state can 
even exceed the one above T c . Obviously, the new triplet 
component generated by the half-metal should provide 
more favorable conditions for the investigation of this un¬ 
usual phenomenon as well as for the observation of the 
related Fulde-Ferrcll-Larkin-Ovchinnikov instabilities— 
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Appendix A: Effective boundary conditions at 
F/HM interface 

Let us derive the effective boundary condition for the 
f tx component of the anomalous Green function at the 
F/HM interface of the S/F/HM heterostructure. From 
Eq. HI) using the condition d x ft y = 0 at x = — ds and 
x = dF + dn one finds the solution for the ft y in the S and 
F layers as well as the solution for in the HM-layer: 

S : fty — A cosh [g s (x + 4)], 

F : f ty =B cosh [qj (x — df)\ + Csinh [qf (x — df )], 

HM : f n = Hcosh[q h (x - d f - d h )], 

(Al) 

where qj = ^/2uj n /Dj (the index j = s, /, h corresponds 
to the S, F and HM layer respectively), Dj is the dif¬ 
fusion constant in the j- th layer, and A, B , C and H 
are the integration constants. Introducing the parame¬ 
ter v sf = {a s /(jf) D f /D s and v h j = (crh/vf) \JDf/D h 
and using the boundary conditions at x = 0 and x = d 
we obtain the systems of equations: 

B cosh(qfdf) — C smh(qfdf) = Acosh(q s d s ), 

— B sinh(qfdf) + C cosh(qfdf) = v s jA sinh (q s d s ) , 
ftx(df) + iB = 0, 

—ftx(df) +iB = H cosh ( q h d h ), 

—d x ftx{df) + iqfC = -qfV h fH sinh (q h d h ). 

(A2) 

Excluding the constants A, B, C and H from this system 
we obtain the effective boundary condition (0) for the 
component ft x - 


Appendix B: Calculation of the anomalous Green 
function in the S/F/HM spin valve 

As follows from Eq. Q if the gap potential A is cho¬ 
sen to be real then the components f s and ft y of the 






anomalous Green function are also real while the compo¬ 
nents f tz and f tx are imaginary. Then it is convenient to 
introduce the complex function 

F(x) = f s + ftzcosd + ftx sm9 (Bl) 

and the real function 


R{x) = i ( f tz sin 9 - f tx cos 9 ), 

so that 

fs = R(n 

ftz = i$s(F) cos 9 — iR sin 9, 
ftx = iSs(F) sin 9 + iR cos 9. 


(B2) 


(B3) 


The introduced functions satisfy the systems of equations 
D. 


2 d x F = K + ih) F — A, 
~J d x R = “nR- 


(B4) 


The boundary conditions for F and R straightly follow 
from the corresponding conditions for the / components. 

Since we assumed that d s -C £ s o the solution of 
Eq. (IBdl) inside the superconductor satisfying the bound¬ 
ary conditions at x = —d s can be represented in the form 


A 

F= — 


Fs 


q 2 s (x + d s y 


R = R s 


q 2 s {x + d s y 


(B5) 


The analogous solution in the F layer reads 
F = JFi cosh [q(x — df)} + F 2 sinh [q(x — df)] , 

R = R\ cosh [qf(x — df)} + i ?2 sinh [qf(x — df)] , 


(B6) 


where we have introduced the complex wave-vector q = 
yj2 (u> n + ih) /Df. 

To calculate the unknown amplitudes in the functions 
F and R let us first consider the boundary conditions at 
x = df. Since f s = ftg = 0 at this interface one finds: 

3?(Fi)=0, 3(fi) cos0 = R\ sin0. (B7) 

These two equations can be rewritten as 

i 7 ! = iRi tan 9. (B8) 

The condition for the component f tx reads [see ©] 

S(qF 2 ) sin0 + g/i ?2 cos0 = — g/T [3(i 7 'i) sin0 + R\ cos0]. 

(B9) 

The boundary conditions at x = 0 give the rest four 
equations: 

A 

-b Fs = F\ cosh (qdf) - F 2 sinh (qdf), 


(BIO) 


where 


P = Psqf 


(j s 2 d s ^ 
a/ D s 


n • 


(BH) 


To solve the system of equations (IB8I) - (1B10I) it is conve¬ 
nient to exclude Fs and Rs from (IB 101) and then express 
R 2 and F 2 in terms of R\ using (IB8I) . The result is [see 
Eq. ©] 


R2 (r 2 fi h ) Ri, 


f 2 = 


iRi tan 9 [pcosh(gd/) + q sinh(gd/)] — pA/u] n 


q cosh {qdf)+p sinh(qd f) 
Finally, substituting (IB12I) into (IB9I) we obtain: 
A ^ s fEsinf?cos6> 


Ri = 


u n Q sin 2 9 + T + (T — 2p,h) cos 2 9 ’ 


where 


W = Q 


qcosh(qdf) + psinh(qdf) J 


(B12) 

(B13) 

(B14) 


J q p + qtanh(qdf) \ 
\q f q + ptanh(qdf) J 


(B15) 


The obtained explicit expression for Ri enables straight¬ 
forward calculation of all other amplitudes in the anoma¬ 
lous Green function. In particular, 


Fs 


q [iR± tan0 — ( A/ui n ) cosh(gd/)] 
q cosh (qd f)+p sinh(q , d f ) 


(B16) 


Appendix C: Calculation of the current-phase 
relation for the S/F/HM/F/S junction 


Let us denote the coordinates of the left and right 
boundaries of the HM-layer as it = — d*/2 and xr = 
dh/ 2. It is convenient to represent the solution of the 
Usadel equation in the Fi-layer in the form 

fs + ftx = Af sinh [q (x - x L )\ + B x cosh [q (x - x L )\, 


fs - ftx = A x sinh [q* (x - x L )\ - B x cosh [q* {x-x L ) 
sinh [qf (x - Xl + df)] 


f ty = iBi- 


sinh (q f d f ) 


ftz = 0, 


(Cl) 

where qf = \J2uj n /Df and q = (co n + ih)/Df (we as¬ 
sume that the diffusion constants and normal conductiv¬ 
ities in the F-layers are equal to each other). In Eq. (1C 1 1) 
we took into account that f s = fn = 0 at x = Xl and 


pFs = —qFi sinh (qdf) + qF 2 cosh (qdf), 

Rs = R\ cosh (qfdf) — R 2 sinh(qfdf), 
pRs = -<l}Ri sinh(cjr/d/) + g/I ?2 cosh(g/d/), 
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ft y = 0 at x = Xl — df. Analogously, the solution in the 
F 2 -layer can be written as 

fa + fty = aJ sinh [q (x - x R )] + B 2 cosh [q (x - x fl )] , 
fa - fty = A 2 sinh [q* (x - £«)] - B 2 cosh [q* ( x - !«)], 
sinh [qf (x-x R -d f )\ 


x = —df — dh/2. Solving the system of equations 
(1C 1 111 we find: 


Pi = 


1 


ie~ i{p / 2 + e i<p f 2 


2 qf coth (qfdf) + 3? [q coth(<?G?/)] 
and, as a consequence, 


sinh(qd/) 

(C12) 


sinh(g/d/) 


ftx = iB 2 
ftz = 0 . 

(C2) 

Inside the HM-layer the solution for the only non-zero 
component /-jq = — ftx + ifty has the form 


*(f:d x fs) = 


i 


sin (ip + 7t/2) 3 2 


sinh(g<2/) 


2 qf coth (qfdf) + 3? [q coth(gd/)] 


(CIS) 


fn = P i sinh (q h x) + P 2 cosh (q h x), 


(C3) 


with qh = u n /D h . Taking derivation of Eq. (IC3I) and 
excluding the constants P± and P 2 we obtain two equa¬ 
tions which connect the values of fff and d x pf on the 
left and right sides of the half-metal (the corresponding 
values are indicated by the upper indexes L and R ): 


Qhfj\ = d x ff t smh(q h d h ) + q h fffCOsh(q h d h ), 
&xf |q — Iff cosh (qitdh ) -f- qhj'ff sinh(('// i d/, ). 


(C4) 


Further for simplicity we will assume that dh <C 
\jDyjA 7rTKo- Then the system (IC4II transforms into 
= f£ j. and d x fQ = d x f^.. Taking this into account 
and substituting Eqs. ca-jca into the boundary con¬ 
ditions at the S/F and F/HM interfaces we obtain: 


Substituting this expression into Eq. m we obtain the 
desired current-phase relation. 


Appendix D: Calculation of the Green function for 
the S/F/HM system of atomic thickness 

Let us introduce the following Green functions in the 
imaginary time representation: 

Goi ,/3 = — (T t (tp a , tp'p'j ^ — (Tt (PaiPp'j ^ J 

[t t ty a , <p+) ) , F++ = (T r ( v £, tf) } , 


E v a = - 


K,P = ~ (Tr (va, tf) ) , ^ (r?+ <fi+) ) 


(Dl) 

Then performing the Fourier transform we obtain the 
following system of the matrix Gor’kov equations: 


— Af s'mh(qdf) + B± cosh (qdf) = — e * V//2 , 

(C5) 

(iu — £)G — tE^ + AIF+ = i, 

(D2) 

— sinh (q*df) — B 1 cosh (q*df) = —e - *^ 2 , 

Mn 

(C6) 

(iu + £) F+ + tP v ’ + - A*IG = 0, 

(D3) 

A 2 sinh(qdf) + B 2 cosh (qdf) = 

n 

(C7) 

(iu -C^jE^—tG- tE v = 0, 

(D4) 

A .2 sinh {q* df) — B 2 cosh (q*df) = —e I¥3 ^ 2 , 

(C8) 

(iu + (?) F^+ + tF+ + tF v+ = 0, 

(D5) 

- B 1 = iB 2 , 

(C9) 

(iu - p) F 7 ' - tE^ = 0, 

(D6) 

- Bi = iB 2 

(CIO) 



—qAf + q* A l — 2Biqf coth (qfdf) = 

(Cll) 

— 2iB 2 qf coth (qfdf) + iqA 2 — iq* A 2 . 

To calculate the Josephson current through the junc¬ 
tion we need to find only the combination (f*d x f s ) at 

(iu + p) F ’»+ + tF^+ = 0, (D7) 

where / = ia y , and 1 is the unit matrix in the spin space. 
Solving the system (ID2MD7I) and considering only the 
linear term in the expansion of the function F + over |A| 
we obtain the expression (1221). 
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